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We consider N fermions in a two-dimensional harmonic oscillator potential interacting with a 
very short-range repulsive pair-wise potential. The ground-state energy of this system is obtained 
by performing a Thomas-Fermi as well as a self-consistent Hartree-Fock calculation. The two results 
are shown to agree even for a small number of particles. We next use the finite temperature Thomas- 
Fermi method to demonstrate that in the local density approximation, these interacting fermions 
are equivalent to a system of noninteracting particles obeying the Haldane-Wu fractional exclusion 
statistics. It is also shown that mapping onto to a system of N noninteracting quasi-particles enables 
D ■ us to predict the energies of the excited states of the N-body system. 
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,__, ■ I. INTRODUCTION 

Haldane |Q first proposed the generalised exclusion 
statistics in the context of excitations in antiferromag- 
netic spin-chains. In this generalised scheme, one consid- 
ers a finite dimensional single-particle Hilbcrt space with 
N identical particles. Increasing the particle number by 
one results in the blocking off of a single-particle states 
for the others, the special cases a — I and a — corre- 
sponding to fermions and bosons respectively. It is to be 
noted that the statistical parameter a is TV- independent, 
and is determined solely by the interaction strength be- 
tween the particles. An exact realization of this statis- 
tics is found in the Calogero-Suthcrland-Moser (CSM) 
one-dimensional model 0, in which the pair interaction 
potential falls off as the inverse square of the relative 
straight-line distance between the particles, which are ei- 
ther on an infinite line, or on a circle. Both the ground- 
state as well as the excited state energies of this interact- 
ing CSM-system may be exactly reproduced by a set of 
N quasi-particles obeying Haldane statistics H. The fi- 
nite temperature distribution function for the generalised 
exclusion statistics is known 0. For the CSM system, it 
has been shown || that the noninteracting quasiparticles 
have such occupancies. In principle, generalised exclusion 
statistics may be obeyed in higher dimensions also. 

The aim of this paper is to present a many-body Hamil- 
tonian in two dimensions that obeys fractional exclusion 
statistics. This model was first proposed in ||, and con- 
sisted of N fermions in a two-dimensional confinement 
potential, and interacting with a very short-range repul- 
sive two-body force. For the one-dimensional CSM case, 
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the ground and the excited state energies of the many- 
body system are analytically known, so one may test the 
mapping to the quasi-particle spectrum. The difficulty 
in the two-dimensional case is that the energy spectrum 
of the proposed many-body system is not known - in 
fact not even the ground state energy is known, so the 
model cannot be tested unless elaborate numerical calcu- 
lations are performed. In ref |j] , only the Thomas- Fermi 
(TF) expression for the spatial density had suggested 
that the system may obey fractional exclusion statistics. 
No attempt was made in that paper to construct a quasi- 
particle scheme to map the energy spectrum. In this pa- 
per, we perform a self-consistent Hartree-Fock (HF) cal- 
culation to estimate the gound-state energy of the two- 
dimensional system. It is found that the simple TF cal- 
culation for the energy agrees surprisingly well with HF 
result. We then use the finite temperature TF method to 
show that in the local-density approximation, the inter- 
acting fermions are equivalent to a system of nonintcr- 
acting particles obeying the Haldane-Wu fractional exclu- 
sion statistics. Finally, we present a quasi-particle spec- 
trum that enables us to calculate the ground as well as 
the excited state energies of the system. Although not ex- 
pected to be exact like the CSM in the one-dimensional 
case, the quasi-particle method yields analogous finite- 
N corrections that lowers the energy. Our model may 
be useful for the description of interacting electrons in a 
quantum dot, since it shows that the short-range com- 
ponent of the effective interaction in a two-dimensional 
system may induce fractional exclusion statistics. The 
mapping to a system of non-interacting quasi-particles 
further enables us to easily calculate the excited states 
of a system of electrons in a quantum dot, albeit ap- 
proximately. This, in the basis of electrons, would be a 
daunting task indeed. 



II. THE MANY-BODY PROBLEM IN 
TWO-DIMENSIONS 

A. The Hamiltonian 

We consider a system of N fermions in two dimen- 
sions, confined in a one-body potential V\, and interact- 
ing pair-wise through a short-range repulsive two-body 
interaction V2. The Hamiltonian is given by 

N N 

F =^E^+E fiM + E ^(Irj-rfcl). (1) 

i— 1 2—1 j<k 

As in ||, the spins are taken to be unpolarized. The 
mean-field expression for the energy at zero temperature 
is given by 



E{p(r)}=Jd 2 



-{v) + V 1 {r)p{v)+ l -{p{v)J d 2 r'p(r')V 2 (\r - r'|) - I / ,/V|,,[v. r' )|-'r j( |r - r'| )} 

(2) 



where r(r), p(r) are the spatial single-particle kinetic 
energy and number density, and p(r, r') is the density- 
matrix. As shown in la], the density-matrix expan- 
sion following Skyrme |j], together with angle-averaging, 
yields the form 



E[p(r)] = Jd 2 



-(r) + li(r)p(r) + ip 2 (r)Af + i{V(r) - [Vp(r)] 2 }M 2 

(3) 



where M n = J d 2 r r n V2{r) is the n th moment of the 
two-body potential. If V^(r) has short range, the higher 
moments fall off rapidly. In this paper we only retain the 
lowest moment Mo, and drop the terms involving M2, M$ 
etc. This is tantamount to taking an effective interaction 
of zero range. In this approximation, we show that the 
interacting system may be mapped onto N noninteract- 
ing quasi-particles obeying fractional exclusion statistics. 
We may therefore consider an effective Hamiltonian 

N N 2 
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(4) 



where a = (1 + ™Tg ) is dimensionless, and will play the 
role of the statistical parameter. Note from the above 
equation that a = 1 describes a system of noninteracting 
fermions. 



B. Thomas-Fermi (TF) approximation 

Consider, within the TF formalism, the above system 
at zero and finite temperatures. We show that this sys- 
tem of interacting fermions is equivalent to a set of nonin- 
teracting particles in V\ (r) obeying fractional exclusion 
statistics locally. The usage of the term locally will be 
explained shortly. Let us first recapitulate the T = 
case. In the TF approximation, r(r) in Eq.ra) is replaced 
by h irp 2 (r)/2m. For the zero-range interaction for V2 
given by (Q), only the lowest moment Mo in Eq.(|3|) is 
nonzero. For this case, (with a spin-degeneracy factor of 
2 included) , the TF energy density functional is given by 



E[p(r)] = Jd 2 r 



fl 2 7T?? 2 

— 7rp 2 (r) + ^(r)p(r) + — (a - l)p(r) 
zm Zm 

(5) 



The ground state energy Eq and the spatial den- 
sity Po F (r) are now determined by the variation 8{E — 
p(0)N) = 0, where /i(0) is the chemical potential at 
T = 0, such that J d 2 r Po F ( r ) — N. On perfoming 
the variation, we get 



m 



Po*(*) = -zs-\m(P)-Vi(t)], r<r , (6) 

irn a 

= 0,r>r , (7) 



where ro is the classical turning point, defined by p(0) = 
Vi(r ). The form of Eq.M for the spatial density 
prompted us to suggest in S that it may also be gen- 
erated by a set of nonintcracting particles, but with a 
zero-temperature occupancy factor 1/a instead of 1. In- 
deed, the same feature was earlier found for CSM, 
where the ground-state energy was also reproduced cor- 
rectly to leading order in N. The chemical potential itself 
is determined by the requirement that the spatial density, 
integrated over space, should give the particle number N. 
For the special choice V\{r) — mu) 2 r 2 /2 which we will ex- 
plore in some detail, this gives the relation 

/x(0) = huVaN. (8) 

For this case, so far as the chemical potential is con- 
cerned, the energy scale has been stretched by a fac- 
tor \fa, an important clue for the construction of 
the quasi-particle spectrum later. On performing the 
r— integration in (S) with the po{r) given by (Q), and 
making use of the above expression for /i, we get the 
ground state energy 

El F = V/2^3/2^ _ (Q) 

O 

Our considerations here are all based on the TF model. It 
is important to test the accuracy of the model, specially 
for small N, by using a more elaborate method. To this 
end, we present in the next section the self-consistent 
Hartree-Fock results. 

Before doing this, however, we demonstrate that this 
system approximately obeys, within the TF frame-work, 
fractional exclusion statistics even at a finite temperature 
T. A system of noninteracting particles at a finite tem- 
perature obeying fractional exclusion statistics has the 
following distribution function at equilibrium tempera- 
ture T 

n a (e,T)= * , (10) 

w(e,T) + a 

where w is the solution of the equation w a [l + w] 1_Q = 
e ( e -M)/7\ p or t ne special case of a two-dimensional trans- 
lationally invariant system (with constant spatial density 
p and constant single-particle density of states), Wu ||] 
had shown that 

4-2 

H{T) = a^—^ +Tln[l - cxp(-7r^ 2 p / 'mT)] . (11) 
m 

In pi, we applied finite temperature TF for a system 
of fermions in a potential Vi(r), and interacting with a 
zero-range force. In this case it is easily shown that 

\u(T) - Fi(r)l = a * h pijC,T ^ + Tln[l - cxp(-7r?i 2 p(r, T)/mT)] 
m 

(12) 
Here p(r, T) is given by M 



p(r,T) = ^ \n[l + exp{(p(T)-V(r,T))/T}} . (13) 



This is the local version of Eq.(|ll|) above, where p(T) 
and a constant p are replaced by their local values 
[p(T) — Vi(r)] and p(r,T) respectively. So, in the pres- 
ence of a fixed external potential and a small number N 
of electrons, as in a quantum dot, for example, fractional 
exclusion statistics may be obeyed in the local- density ap- 
proximation. The latter term means that at any point r, 
the local fermi-energy and spatial density obey the same 
equations as their counterparts in the translationally in- 
variant system. This generally is true in the interior of 
the system, but deviations may occur at the edges if the 
spatial density changes rapidly in a scale compared to 
the size of the system. The local-density approximation 
has been successfully used in sophisticated nuclear many- 
body calculations even with a relatively small number of 
nucleons ||, where a large fraction of them may reside 
in the surface region. 

In case the single-particle density of states is not a 
constant, the Wu-equation (|ll| ) is generalised to the fol- 
lowing form |Q: 

p{T) = av + Tln[l - exp(-u/T)] , (14) 

where 



i/(T) = I n a (e,T)de. (15) 

Jo 



We emphasize that Eq.(14) is quite general, and does 



not make any restrictive assumption about the density 
of states. We now show that our interacting fermionic 
system also obeys a local version of Eq.pl) in the TF 
frame-work. To demonstrate this, we define 

u(t,T)= / den F {r,p,T) , (16) 

Jo 

where 

n F (r, p, T) = -^— — wtTTI ' ( 17 ) 

exp[(p z /2m + V(r) — p)/l \ + 1 

For the system given by the effective Hamiltonian (0), 

V(r,T) = V 1 (r) + —(a~l)p(r,T). (18) 

m 

Replacing de in Eq.(Jlq) by p dp/m at a fixed point r, and 
performing the p-integration, we then obtain 

„(r,T) = ^_p(r,T). (19) 

m 

Using Eqs ( jTg) ) and (12), both based on the interacting 
fcrmion system, we can then write 

[p(T) - Vi(r)] = ai/(r, T) + Tln[l - exp(-v(r, T)/T)\ . 

(20) 



This is the local version of Eq.(|Tj) that we set out to 
show, and emphasizes that the local density approxima- 
tion is applicable even when the single-particle density of 
states is not a constant. 

We now go on to check the accuracy of the TF ap- 
proximation for the ground state by performing a self- 
consistent HF calculation in the next section. 



C. The Hartree-Fock (HF) calculation 

The advantage of performing a HF calculation with 
H given by (Q) is that the Fock term is local. A little 
algebra shows that the HF single-particle orbitals ipi(r) 
with energies et obey the one-body equation 
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h + (a-l)p(r) 

m 



^i(r) = eiV>i(r) , (21) 



where ho = p 2 /2m + mu> r /2, and p = J^ V'*V'j t ne 
prime indicating that the sum is over the occupied or- 
bitals only. Note that the one-body potential in (£l|) is 
the confining harmonic oscillator in ho, plus the contri- 
bution coming from the direct and the exchange matrix- 
elements of the two-body interaction. Let us denote the 
latter by 

[/(r) = — (a-l)p(r). (22) 

m 

The potential U, proportional to p, is determined by the 
occupied single-particle orbitals. Therefore the equations 
have to be solved iteratively, starting with a trial p, until 
self-consistency is achieved to a specified accuracy. To 
this end, the orbitals ipi (r) are expressed in the basis set 
generated by the unpcrtubcd harmonic oscillator. The 
basis set, of course, has to be truncated in practice. For 
numerical accuracy, we find that if M shells are occu- 
pied in the unperturbed configuration, then at least 2A/" 
shells should be included in the basis set. For N = 12 
particles, where J\f = 3, we have tested that the self- 
consistent results remain essentially unaltered when the 
dimension of the basis set is increased from shells 6 to 
10. The HF single-particle Hamiltonian, given by the 
left-hand side of Eq. (f2l|) (with a trial p), is represented 
by a matrix in this basis. It is diagonalsied to obtain 
the orbitals and their eigenvalues. The density p and 
the corresponding HF Hamiltonian are recalculated, and 
the procedure repeated for self-consistency Typically, for 
moderate strength of the interaction (a = 3/2), about 40 
iterations or more are needed to obtain self-consistent oc- 
cupied orbitals tpi- The HF ground-state energy is given 
by 



Er=Y,U-\<^mi>) - (->:-!! 



where U is defined by Eq. (E2J) , and the sum over i is over 
occupied states only. The self-consistent results for 3 dif- 
ferent systems (N = 12, 20 and 42 fermions), each for 
2 values of a, are presented in Table 1. These N- values 
are closed shells configurations for ho. For non-closed 
shell systems, we found that numerical convergence for 
self-consistency is harder to achieve. Interestingly, as 
is shown in Table 1, Eq F is very close to E^ F calcu- 
lated from Eq.(H). In Fig.l we display, for N = 20 and 
a = 3/2, the shape of the self-consistent HF one-body 
potential, {muj 2 r 2 /2 + U(r)) as a function of the dimen- 
sionless variable R — (muj/fiY^r, and the corresponding 
low-lying single-particle energies e^ in units of Tiuj. For 
comparison, the unperturbed harmonic oscillator and its 
energy eigenvalues are also shown in the same diagram. 
Two features are to be noted : (a) the relatively large 
upward shift at the center in the HF one-body potential, 
resulting in an overall shift in the single-particle spec- 
trum, and (b) the almost equal but reduced spacing be- 
tween the levels in the HF spectrum, with the degeneracy 
between the (2s, Id) states being intact to 4 significant 
figures, and the splitting in (2p, If) and (3s, 2d, lg) sub- 
shells remaining negligble. These characteristics suggest 
that, like the CSM, the ground and the excited states 
of the N-particle system are shifted by the same overall 
constant whose magnitude is determined by the stength 
of the interaction. (For systems in which the last filled 
shell is only partially occupied, preliminary results indi- 
cate that the splittings between the states within a shell 
are larger.) We also note that the states (n, I) and (n, —I) 
remain strictly degenerate. We also found that in every 
case, the spacing between the HF shells are almost equal, 
with an effective oscillaor constant given by ui/y/a. This 
behaviour of the HF e^ 's is very different from the quasi- 
particle spectrum e;'s presented in the next section, the 
latter actually spreading out by a factor y/a. We may 
express the HF single-particle energies in the form 

e„ = A(a,Af ) + (n+l)u>/Va ,n = 0,l,.. (24) 

where A/o is number of filled shells in the system. The 
gap A is proportional to (A/o + 1), and the constant of 
proportionality increases with a. 

As mentioned in the previous section, the TF result 
(fjj) had suggested that the ground state energy of the 
interacting system may be estimated by filling up each 
level of the corresponding noninteracting system by a 
fraction 1/a, but preserving the total particle number. 
We call this the "global" fractional-statistics rule. It is 
useful only for the ground state, in contrast to the quasi- 
particle formalism given in the next section. In Table 
1, the numerical values for E G (global) are also given for 
comparison. 



III. THE QUASI-PARTICLE SPECTRUM 

We wish to express the energy of an interacting N- 
particle system as a sum of noninteracting quasi-particlc 



energies, both for the ground and the excited states. This 
can be done exactly for CSM, which behaves as an ideal 
gas obeying fractional exclusion statistics. For the two- 
dmensional model presented here, we cannot prove that 
such is the case. Nevertheless, we show that when the 
degeneracies of the levels are taken into account, a con- 
sistent scheme for the quasi-particle spectrum may be for- 
mulated, obeying the same rules as given for CSM ||U . 
This procedure reproduces the large-N TF results for the 
ground state, and predicts the finite-N corrections. Fol- 
lowing the steps of 0], we then show that such a map- 
ping leads naturally to the thermodynamic distribution 
function formulated in B. 

In what follows, we take the two-dimensional har- 
monic oscillator as the confining potential, and a spin- 
degenaracy of two. At the outset, we clarify certain def- 
initions that are used in this presentation. The unper- 
turbed harmonic oscillator spectrum is specified by shells, 
containing degenerate states in the shell. For example, 
the j th shell is at energy ju> and has r)j = 2j states in it, 
where we have included the spin degeneracy of two in rjj . 
The quasi-particle spectrum has exactly the same struc- 
ture of these shells and states, except that the energy of 
a quasi-particle is determined by the occupancies of the 
states below it, as specified by the rules below. In this 
mapping, the unperturbed energy e k = koj of a particle 
is mapped onto the quasi-particle energy e^(fc,a). The 
analogy with CSM leads us to suggest the following 

e A (k,a)=e k -uj(l-V^W{e k ,0), (25) 

where Af(e k , 0) is the sum of the occupancy fractions of 
the shells with energies e < efc. This is given by 

M{t k ) = J2 9(e fc - ej)- J2 n jti , (26) 

j=i "* i=i 

where Q(x) = for x < and 1 for x > 0, and rijj 
is the occupancy (0 or 1) of the i th state of the j th 
shell. An example is illustrated in Fig. 2, where (a) is 
an unperturbed harmonic oscillator excited state con- 
figuration in which one particle from the (2s, Id) shell 
has been excited to the next (2p, 1/) shell. In (b), the 
corresponding mapping to the quasi-particle spectrum is 



shown. We use the notation Q — y/auj. Using Eq.(|26[), 
we note that A/"(ei) = 0, A/(e 2 ) = 1, A/"(e 3 ) = 2jy[e 4 ) = 
17/6, Af(e 5 ) = 71/24. Now we may use Eq.(§|) to de- 
duce the quasi-particle energies as shown in Fig. 2(b). 

For the ground state of a closed-shell configuration, all 
the shells are fully filled upto (say) the A/'q' 1 shell, and 
are unoccupied thereafter. The above rules then lead to 

Af(e j )=j-1, l<j<A/" ; 

= A" , j>A/" . (27) 

The energy levels e^(A:, a) in the ground-state configura- 
tion are given by ui,u> + uj,u> + 2cj, ...u> + (A/'q — l)w, 2w + 

(A/o — 1)&, 3uj + (A/"o — l)w, The ground state energy 

is then given by 



4 9P) =f>(fc,a) % , 
fe=i 

A/"o jVo 

= 2^fce fe -2w(l- Va) 52 A(fe - 1) , 

fe=i fe=i 

= ^V^AA (AA + l)(2A/b + l)/3 - w(V^ - 1)AA (AA + 1) . (28) 

where Ao is the last filled level. Expressing Ao in terms 
of N, 

N = 2j2k, (29) 

fe=i 



the ground state energy is given by 



E^ qp) = lui^N 3 ' 2 



1 

8iV 



~(y/a-l)uN . 
(30) 



The term proportional to TV 3 / 2 is the same as the TF 
expression ([)]), while the others are the finite-N correc- 
tions arising from the microscopic approach taken here. 
The values obtained for E^' by Eq.(pOl) are displayed 
in Table 1. Note that these are always lower than the 
corresponding HF values, as is to be expected. Note also 
that for a = 0, Eq.(|30|) gives the bosonic limit with all 
particles in the lowest energy state at ui. Although the 
expression (|3C]) is only valid for a closed-shell system, the 
ground state energy of an "open shell" configuration may 
also be calculated easily using this formalism. 

The excited state energies for any configuration {nk,i} 
may now be written using Eqs.(E5l f26[): 



E[{n k ,i}} = *^2 e A(k,a) 52 n M i ( 31 ) 

fe=l i=l 

oo Vk oo oo , Vk Vk 

= ^2 e k^2n k s-uj {I- yfa) 52 52 — 52 '^2 n k,t n k',i' ■ (32) 
fe=l i=l k'(<k) = lk=l ^ k ' i=l i' = l 

Note that the energy levels e^(fc, a) depend on the con- 
figuration {nk,i}- Like the CSM, the excitation energies 



of the system are given only by the first term in Eq.(32). 
The second term depends only on the total number of 
particles. The correspondance with CSM and the agree- 
ment of the leading N-dependent term of Eq P with the 
TF and the HF results suggests that (|3^) may give a good 
description for the excited states. As an aside we remark 
that if the effective interaction between N electrons in a 
quantum dot is dominated by a short range interaction, 
then the above expression for the total energy as a sum 
of the quasiparticle energies immediately yields not only 
the approximate ground state but also its excited states. 
A direct calculation of these excited states is in general 
extremely difficult. 

We now consider the distribution function for the 
quasi-particles. First consider the ground state (T = 0). 



So far as the energy spacing between the occupied quasi- 
particle levels are concerned, it is given by Co = y/au/. 
The number of particles ANa in an energy interval be- 
tween €a and €a + Aca is given by 

AN A = ~-e A Ae A , 

= kcaAca ■ (33) 

a u) z 

Thus the occupancy factor for the filled quasi-particlc 
states in the energy scale of u> is 1/a instead of 1. This 
is precisely the T = distribution function for the occu- 
pied states in fractional exclusion statistics, irrespective 
of dimensionality. 

In conclusion, we have shown that for fcrmions in a 
two-dimensional harmonic potential interacting with a 
very short-range repulsive two-body force, the TF calcu- 
lation for the ground state energy agrees well with the 
corresponding HF result. The HF self-consistent single- 
particle states get almost a constant upward shift, while 
still retaining equal spacings in a slightly shallower po- 
tential. Gaining confidence with this success of the TF 
method, we apply it for the finite temperature problem 
in a confinement potential V\(r). We find that in this 
approximation, the interacting fcrmions obey the same 
equations locally as to be expected from a system of non- 
interacting particles obeying fractional exclusion statis- 
tics. This is true for any confinement potential, so long as 
the finite-range effects of the two-body interaction is ne- 
glected. It means that for such a system, the interacting 
fcrmions in the bulk obey fractional exclusion statistics, 
atleast so far as energies are concerned. For the special 
case of harmonic confinement, a quasi-particle spectrum 
is constructed to predict the ground- as well as the ex- 
cited state energies of the system. This is analogous to 
the CSM case in the one-dimensional problem. 
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TABLE I. Ground state energy of N-particle (closed shell) 
systems for various iV's (column 1) and interaction strengths 
(column 2) . Column 3 displays the sum of the occupied or- 
bital energies, and the next the HF ground state energy (see 
Eq.(H). Column 5 is the TF energy given by Eq.@. Col- 
umn 6, with the heading E , is the "global" approximation 
as explained at the end of the HF section. The last column 
is calculated from the "microscopic" quasi-particle approach, 
as given by Eq.feOJ). 



N a %ei E HF E TF E a B» 



12 


5/4 


34.424 


31.301 


30.984 


32.000 


29.890 


12 


3/2 


39.969 


34.278 


33.941 


34.667 


31.598 


20 


5/4 


73.781 


67.079 


66.667 


68.000 


64.723 


20 


3/2 


85.694 


73.472 


73.030 


73.333 


68.991 


42 


5/4 


223.818 


203.478 


202.879 
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FIG. 1. Plot of the HF self-consistent potential 
(mu) 2 r 2 /2 + U(r)) and the corresponding single-particle en- 
ergies e; (in bold), for N — 20 and a = 3/2. For compari- 
son, Vi(r) — mui 2 r 2 /2 and its energy levels (in light) are also 
shown. The distance R is in units of ^Jh/mui, and the energy 
scale is in units of hui. 

FIG. 2. The quasi-particle spectrum for a one-particle 
one-hole excited state configuration for N = 12. The crosses 
denote the particles, and the circle the hole. On the left, 
the configuration is shown for the harmonic potential (unper- 
turbed), to the right is the quasi-particle spectrum to which 
the interacting system has been mapped. 
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